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Abstract 

We calculate differential cross sections for exclusive double diffractive (EDD) production of 
open charm in proton-proton and proton-antiproton collisions. Sizeable cross sections are found. 
The EDD contribution constitutes about 1 % of the total inclusive cross section for open charm 
production. A few differential distributions are shown and discussed. The EDD contribution falls 
faster both with transverse momentum of the c quark/antiquark and the cc invariant mass than 
in the inclusive case. 
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I. INTRODUCTION 



The open charm production is often considered as a flag reaction to test the gluon dis- 
tributions in the nucleon. For the cc and bb production at high-energies the gluon-gluon 
fusion is assumed to be the dominant mechanism. This process was calculated in the NLO 
collinear [l| as well as in the fc^-factorisation 0, [s], 0, @| approaches by several authors. These 
analyses seem to report on missing strength^. This suggests that other processes ignored so 
far should be carefully evaluated. 

The number of potential contributions is not small. In the present paper we concentrate 
on exclusive double diffractive (EDD) mechanism, which was not considered so far for the 
cc production. The mechanism of the exclusive double-diffractive production of open charm 
is shown in Fig. [H 




FIG. 1: The mechanism of exclusive double-diffractive production of open charm. 



The EDD bb reaction consitutes a irreducible background to the exclusive Higgs boson 
production [gI] measured in the bb channel. Up to now only approximate estimates of the bb 
production were presented in the literature. In the present paper we consider the pp ppcc 
reaction as a genuine 4-body process with exact kinematics which can be easily used with 
kinematical cuts. The amplitude of the genuine four-body reaction is written in analogy 
to the Kaidalov-Khoze-Martin-Ryskin (KKMR) approach used previously for the exclusive 
Higgs boson production 0, @, 0] • 



II. MATRIX ELEMENT AND THE CROSS SECTION 

FOR EXCLUSIVE DOUBLE DIFFRACTIVE qq PAIR PRODUCTION 

Inclusive heavy quark/antiquark pair production was considered in detail, e.g. in Refs. 0, 
The nonrelativistic QCD methods were successfully applied also in the case of central 
exclusive production of heavy quarkonia in Refs. [l^, ll|, 12]. It looks quite natural to apply 
similar ideas to exclusive diffractive qq (unbound) pair production. 



The situation is often somewhat clouded by studying the uncertainty bands due to variation of renor- 
mahzation and factorization scales. These analyses lead to rather broad uncertainty bands which prevent 
definite conclusions. 



2 



A. Kinematics 



The kinematical variables for the process pp p + "gap" + 
in Fig. H 



+ "gap" + p are shown 
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FIG. 2: Kinematical variables of exclusive diffractive production of qq pair. 

We adopt here the following standard definition of the light cone coordinates 

k+ = n^k"" = k^ + k\ k- = n-k" = k^ - k\ kt = (0, k\ 0) = (0, k,0), 

where are the light-cone basis vectors. In the c.m.s. frame 

P2 _ Pi 
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(2.1) 



and the momenta of the scattering hadrons are given by 

Vt =P2 = Pi =Pt = Pl,t = P2,t = 0, 

with the Mandelstam variable s = AE^^^. 

Within the standard fcrfactorisation approach, the decomposition of gluon momenta into 
longitudinal and transverse parts in the high-energy limit is 



qi = xipi + qi^t, q2 = X2P2 + q2,t, < Xi,2 < 1, 
qo = x[pi + X2P2 + qo,t, ~ 4 < a;i,2, qo,i,2 - qo/i/2,f 

Making use of energy-momentum conservation laws 

qi=Pi- p'l - go, q2 = P2- P2 + qo, qi + q2 = ki + h 

we write 



(2.2) 



2 , l-D 12 _ ,^2 Ml^={ki + k2f 



(2.3) 



(2.4) 



where Mqq is the invariant mass of the qq pair, and is its transverse 3-momentum. 
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B. The amplitude for pp — > ppQQ 

Let us concentrate on the simplest case of production of qq pair in the color singlet state. 



Color octet state would demand an emission of an extra gluon [13| which considerably 
complicates the calculations, and we postpone such an analysis for future studies. 

In analogy to the Kaidalov-Khoze-Martin-Ryskin approach (KKMR) 0, H, |9| for Higgs 
boson production, we write the amplitude of the exclusive diffractive qq pair production 
pp p{qq)p in the color singlet state as 



Qlt lit (lit 



(2.5) 



where Ag, Xq are helicities of heavy q and g, respectively. Above ff^ and are the off- 
diagonal unintegrated gluon distributions in nucleon 1 and 2, respectively. They will be 
discussed in a separate subsection below. 

The longitudinal momentum fractions of active gluons are calculated based on kinematical 
variables of outgoing quark and antiquark 

xx = ^^exp(+?/3) + ^^exp(+?/4) , 

X2 = ^^exp(-?/3) + ^^exp(-?/4) , (2.6) 

where m^^t and rrii^t are transverse masses of the quark and antiquark, respectively, and 7/3 
and 1/4 are corresponding rapidities. 

The bare amplitude above is subjected to absorption corrections which, in general, depend 



on collision energy and on the spin-parity of the produced central system We shall 

discuss this issue shortly when presenting our results. 



C 99 ^ QQ vertex 

Let us consider the subprocess amplitude for the qq pair production via off-shell gluon- 
gluon fusion. The vertex factor V^^^^_ = V^^^^_{qi,q2,ki,k2) in expression (12. 5p is the pro- 
duction amplitude of a pair of massive quark q and antiquark q with helicities Xq, Xq and 
momenta ki, k2, respectively. Within the QMRK approach p^], the color singlet qq pair 
production amplitude can be written as 

V-;-(gi, q2, k,, k2) ^ n;n-V^X'''\q,, g^, k,, A;^), (2.7) 
V^X^'^iqu q2, fci, k2) = -g' J2 (3^' >< 

i,k 

^x,{kim]t%h^''{qi, 92, A^i, k2) - tDf^fb^^iq^, q^, ^i, k2))v^,{k2), 

where are the color group generators in the fundamental representation, u{ki) and f (^2) 
are on-shell quark and antiquark spinors, respectively, 6, h are vertices (12.81) arising from the 
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Feynman rules : 

b^^iqi, q2, h, h) = Yj^^^^^f^Y , (2.8) 

[Qi — 1^2) — 

The SU(3) Clebsch-Gordan coefficient (3z,3A;|l) = 5*7 in Eq. ([22D projects out the 
color quantum numbers of the qq pair onto the color singlet state. Factor l/y/TT^ provides 
the averaging of the matrix element squared over intermediate color states of quarks. 

The tensorial part of the amplitude is therefore: 

.rt^u I , , ^ 2- n ^( u qi-ki-m ^ qi-k + m ^\ 

V,^,^{q„q2,h,k2)=g,u,^{h)^^ ^^^_j,^y_^2 ^'-l \^^_j,^y_^2 ^ J-A.(^2). 

(2.9) 

Taking into account definitions (12. ip and momentum conservation (12.31) and using the 
gauge invariance properties we get the following projection to the light cone vectors (so 
called "Gribov's trick") 

yxX^<li,l2,ki,k2) = n+n-V^'^^^l^^{qi,q2,ki,k2) 



^V..^Wi'^2,A^i,/c2) 



S Xi X2 

Using (12. 9p and (I2.10p we can write 



' 'nX.^(^i'^2,A:i,A;2). (2.10) 



J J ^ 4 2- n \ f- qi-ki-m ^ ^ qi - h + m \ 

(2.11) 

The coupling constants 9s{fJ'r,i)9s{fJ'r,2)- In the present calculation we take the renor- 

mahzation scale to be fi^^ = /i^2 = M^^/A. The matrix element (12.111) is then calculated 
numerically. Inserting it to Eq. (12.51) we can calculate numerically the whole amplitude for 
the pp — > ppQQ process. 
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D. OfF-diagonal unintegrated gluon distributions 



In the KMR approach the off-diagonal parton distributions are calculated as 

.KMR. ^2 „2 , N p d[g{xuk^)Si/2{klfi^)] 

Ji K^uQiplJ^ = Rg dhgkf lfc?=Q?t^l^i) 

Fit,), 



.KMR. ^2 , ,2 , ^ p d[9{x2, k^) Si/2{kl fl^)] 
h {X2,(^2,uf^ yh) - rig ^ ^2 \k^,=Ql^r[t 



2) 



(2.12) 



where Si/2{qf, n ) is a Sudakov-like form factor relevant for the case under consideration [Ti 
The last approximate(!) equalities come from the fact that in the region under consideration 
the Sudakov-like form factors are somewhat slower functions of transverse momenta than 
the coUinear gluon distributions. While reasonable for an estimate of gluon distribution it 
may be not sufficient for precise calculation of the cross section. It is reasonable to take a 
running (factorization) scale M2-/4. We shall call the formulae fl^TT^ as the 

DDT-like formulae [18|, for brevity. 

The factor Rg here cannot be calculated from first principles in the most general case of 
off-diagonal UGDFs. It can be estimated in the case of off-diagonal coUinear PDFs when 
x' <^ X and xg = — a;)" j20|. Then 



"^^-^ r(A + 4) • ^2.13) 

In a more realistic case of DGLAP GDF A = A(x, /i^). Typically i?g ~ 1.3 - 1.4 at Tevatron 
ener gy. A more general case of unintegrated off-diagonal distributions was discussed in 
Ref. 17|, but we will not touch them here. 

The off-diagonal form factors are parametrized here as: 

F(t) = exp (fiofft) . (2.14) 

In practical calculations in this letter we take B^s = 2 GeV~^. 

In the original KMR approach the following prescription for the effective transverse mo- 
mentum is taken: 

Qlt = min(g^^t,g^J , 

Ql, = mm{ql,,ql) . (2.15) 



Other prescriptions are also possible [IC 

In evaluating fi and /2 needed for calculating the amplitude (12. 5p we use the GRV 
coUinear distributions [l6l |. 

It was proposed to express the 5*1/2 form factors in Eq. (I2.12p through the standard 
Sudakov form factors as: 

S,/2{qlfi') = jTg{ql^^). (2.16) 
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FIG. 3: The Sudakov form factor as a function of transverse momentum squared and a logio of 
the scale parameter /x^. 

The Sudakov form factor, a two-dimensional function Tg{q^, function of transverse 

momentum squared qf and a logio of the scale parameter /i^, is shown in Fig. [3l 

A strong dependence on /x^ (log^Q(/i^) in the figure) is clearly visible. This dependence 
leads to a huge perturbative damping of the off-diagonal UGDFs (and, as a consequence, of 
the amplitude and the cross section) in the case when objects X with sizeable masses (Mx) 
are produced, i.e. when fi^ ~ M|^. 



E. The pp ppQQ cross section 

The cross section is obtained by assuming a general 2^4 reaction: 

1 , . , , ^ d^Pi d^P2 d^Ps 



d^P4 



da = -|^,^,P(2.)^5^(;,, -p, -p, -P4)-^2.^32i,, i2.f2E, (2nY2E, ' 

(2.17) 

The details how to conveniently reduce the number of kinematical integration variables are 
given elsewhere 19 . 
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III. RESULTS 



As for the exclusive production of Xc mesons 0,113, in the case when the KMR UGDFs 
are used an extra cut-off on transverse momenta is apphed i.e. the formula (12. 5p is used if 

Q\,t = min(g2^,g2 j > gj^^^ ^ 

Ql,t = min(g^^^,gjj > gj^^^ _ (^^ -^^^ 

Otherwise the cross section is set to zero. 

Let us proceed now with the presentation of differential distributions of charm quarks 
produced in the EDD mechanism. In our calculation here we fix the scale of the Sudakov 
form factor to be /i = Mcc/2. Such a choice of the scale leads to a strong damping of the 
situations with large rapidity gaps between c and c ^. 

In Fig. m we show distribution in rapidity. The results obtained with the KMR method 
are shown together with inclusive gluon-gluon contribution calculated as in Ref. @. The 
effect of absorption leads to a damping of the cross section by an energy-dependent factor. 
For the Tevatron this factor is about 0.1. If the extra factor is taken into account the EDD 
contribution is of the order of 1% of the dominant gluon-gluon fusion contribution. 

The corresponding rapidity-integrated cross sections are: 6.6 /ib for exact DDT formula, 
2.4 /ib for simplified DDT formula (see Eq. (12.121) ). For comparison the inclusive cross 
section (gluon-gluon component only) is 807 /ib. 




FIG. 4: Rapidity distribution of c or c. The upper curve is for inclusive production in the k 



factorization approach with the Kwieciiiski UGDF and ^'^ = 4mc, while the two lower lines are for 



the EDD mechanism for the KMR UGDF with leading-order collinear gluon distribution [16(] . The 
solid line is calculated from the exact DDT- like formula (see Eq. (j2.12p ) and the dashed line for the 
simplified formula (when only derivative of collinear GDF is taken) . An extra cut on the momenta 
in the loop Qt^ut — 0-^6 GeV^ was imposed. Absorption effects were included approximately by 
multiplying the cross section by the gap survival factor 5 = 0.1. 



Large rapidity gap means automatically large invariant mass Mcs of the cc systei 



8 



In Fig. [5] we show the differential cross section in transverse momentum of the charm 
quark. Compared to the inclusive case, the exclusive contribution falls significantly faster 
with transverse momentum than in the inclusive case. 




P,(GeV) 

FIG. 5: Transverse momentum distribution of c or c. The other details are the same as in Fig. [H 

In Fig. E] we show the distribution in the invariant mass of c and c. The fluctuations 
visible in the figure are due to the fact that the integration is not directly in Meg, but in 
other variables, and the number of integration points is rather restricted. Compared to 
the inclusive case the invariant mass distribution for the FDD component is significantly 
steeper. This is due to the Sudakov-like form factor which, according to the procedure 
described above, damps the cross section for large invariant masses M^. 



=' ' ' I ' ' ' I ' ' ' I ' ' ' I ' ' ' I ' ' ' I ' ' ' I ' ' ' I ' ' ' I ' ' ' i 




(GeV) 

FIG. 6: Invariant mass distribution of the cc pair. The other details are the same as in Fig. HI 
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As in the inclusive case within the fci-factorisation approach j5|] the cc pair possesses the 
transverse momentum different from zero. The corresponding distribution is shown in Fig. [71 
The distribution for the exclusive case (the two lower lines) is much narrower compared to 
the inclusive case (the upper line). 




p (GeV) 

t.sum 



FIG. 7: Distribution in the transverse momentum of the cc pair. The other details are the same 
as in Fig. SI 



IV. CONCLUSIONS 

In the present letter we have evaluated, for the first time in the literature, the contribution 
of exclusive-double diffractive production of open charm. We have found a sizeable cross 
sections of the order of 1 % of the standard inclusive gluon-gluon fusion contribution at the 
Tevatron energy. The details depend, however, on the UGDFs used in the evaluation of the 
cross section. The most reliable estimate is obtained with the KMR off-diagonal UGDFs. 
These distributions were verified recently in the production of Xc quarkonia 11, T2|, i-e. for 
the kinematics similar to the present one. 

It would be therefore very valuable to measure the diffractive mechanism discussed in 
the present paper. How to identify the FDD contribution? The events corresponding to 
the FDD contribution are expected to be related to a rather small multiplicity of particles 
(mainly pions or kaons) associated with c or c, or more precisely charmed mesons which are 
formed in the process of hadronization of charm quarks into charmed mesons. One method 
would be therefore to measure D mesons with a trigger on small pion multiplicity. Another 
method would be to measure D mesons in association with rapidity gaps with respect to 
the outgoing protons/antiprotons. This is partially possible at the Tevatron and will be 
accessible at the LHC as well. 
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